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SPH Simulations of Negative (Nodal) Superhumps: A 
Parametric Study 
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ABSTRACT 

Negative superhumps in cataclysmic variable systems result when the accretion disc 
is tilted with respect to the orbital plane. The line of nodes of the tilted disc precesses 
slowly in the retrograde direction, resulting in a photometric signal with a period 
slightly less than the orbital period. We use the method of smoothed particle hydro- 
dynamics to simulate a series of models of differing mass ratio and effective viscosity 
to determine the retrograde precession period and superhump period deficit e_ as a 
function of system mass ratio q. We tabulate our results and present fits to both e_ 
and e+ versus q, as well as compare the numerical results with those compiled from 
the literature of negative superhump observations. One surprising is that while we find 
negative superhumps most clearly in simulations with an accretion stream present, we 
also find evidence for negative superhumps in simulations in which we shut off the 
mass transfer stream completely, indicating that the origin of the photometric signal 
is more complicated than previously believed. 
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1 INTRODUCTION 

Cataclysmic variable (CV) binaries typically contain a white 
dwarf primary Mi accreting matter from a l ow-mass main- 
sequence seco ndary M2 via an accretion disc l|Warner1ll995al ; 
' iHellierl l200lh ■ The Roche-lobe filling secondary star loses 
mass through the LI inner Lagrange point at a rate in the 
approximate range M ~ KT 11 to ~1CT 7 M Q yr -1 . The 
observed system characteristics depend most strongly on the 
mass ratio and individual masses, inclination, mass-transfer 
rate and composition, and the strength and geometry of any 
magnetic fields that may be present. 

Systems with accretion flows that are not substan- 
tially affected by stellar magnetic fields and with mass 
ratios q = M2JM1 in the ra n ge 0.03 <<?< 0.35 (e.g., 
iBaker-Branstetter fc Wood! Il999l ; iMontgornervf l200ll ; this 
work) can have the outer disc extend to radii approximately 
between that of the 3:1 corotation radius and the 2:1 corota- 
tion resonance allowing driving by the m ode-coupling mech- 
anism first identified by iLubowl |l99ll ). Within this range, 
common superhump oscillations can be driven and yield 
photometric time-series with amplitudes of up to ~0.1 mag- 
nitudes and periods slightly in excess of the orbital period. 



The observed fractional period excess of the su- 
perhu mp oscillations i ncrea ses with increasing mass 
ratio (|Patterson et al. I 120051 ). The hydrodynamics of 
the superhump oscillations have been simulated nu- 
merically and extensively with particle-based schemes 
including smoothed particle hydrodynamics (SPH; 
Whitehurstll 19881; iHirose : fc Osakil[l990l; IWhitehurst fc Kind 



19911; IWhitehurstJ Il994l; [Murray! Il99rj, Il998l: 



Simpson fc Wood! Il998l: IWood. Montgomery, fc Simpso: 



2000) ; iFoulkes. Haswell. fc Murray! 120061 ; ISmith et al 

2007) and recently with a 2D grid-based approach 



( Kiev. Papaloizou. fc Ogilvie] 120081 ). First recognized as 
photometric variations with a period a few percent longer 
than the orbital period that occurred during the s uper- 
outb u rsts of the S U UMa subclass of dwarf novae (|Vogtl 
ll974l ; IWarneilll975l ), common superhump oscillations have 
since b een identified in novalike CVs dPatterson et al 



1993bl ; iRetter. Leibowitz. fc o"fek1 Il997l ; ISkillman eta! 
1997), the double helium white dwarf AM CVn bina- 



ries (iPatterson. Halpern. fc Shambrookl Il993al ; IWarnerl 
Il995bl; iNelemand 120051) 7 and in low-mass X-ray b inaries 
JCharles et all Il99ll;~ iMineshige. Hirose, fc Osakil Il992l ; 
lO'Donoghue fc Charleslfl996l ; iHvnes et al.ll2006| ).~ 
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Common superhumps are interchangeably referred to 
as positive or apsidal superhumps in the literature, or 
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just superhumps with no modifier. The term "positive" is 
empirically-based and refers simply to the sign of the posi- 
tive period excess e+ defined as e+ = (P+ — P or b)/ Porb where 
we adopt the notation P+ to represent the positive super- 
hump oscillation period. The term apsidal is dynamically- 
based and refers to simulation results indicating that while 
the radial disc profile as a function of co-rotating azimuthal 
angle is non-stationary, the oscillation-averaged profile is ec- 
centric with a line of apsides that precesses slowly in the 
orbital prograde direction with a period Pp [cc given by 
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(i) 



An accretion disc which is tilted out of the or- 
bital plane will precess in the retrograde direction, 
analogous to the retrograde precession of the line of 
nodes of th e moon's orbit. We first stud i ed th is nu- 
merically in IWood. Montgomery, fe Simpson! (|2000l . here- 
after Paper I, and see also Montgomery 2004). In an 
early paper discussing the eclip s ing b inary TV Col, 
iBonnet-Bidaud. Motch. fc MouchetJ (| 19851 ) suggested that 
the observed 5.2-h period was the difference between the 
orbital frequency ^ or b = l/P or b (where P or b = 5.5 
h) and the natural retrograde precession frequency ^j7 roc 
(fprec ~ 4 d) of a tilted disc. In a follow-up study, 
iBarrett. O'Donoghue. fc Warned (119881 ) suggested that the 
physical source of the observed periodicity was the sweeping 
of the accretion stream bright spot across th e face of a tilted 
disc, a model we were able to confirm in IWood fc Burkel 
(120071. hereafter Paper II) and wh ich was also discussed in 
iFoulkes. Haswell. fc Murray! l|2006l ) . These photometric vari- 
ations - with periods slightly less than P or b and a completely 
different physical origin than the common superhumps - 
have also have come to be classified as superhumps. While 
the name nodal superhumps has been proposed as a phys- 
ically descriptive name for these variations (and is the one 
we prefer), it is most common in the literature to refer to 
them as negative superhumps, from the sign of the relation 
e_ = (P_ — P or b)/Porb, where P_ is the observed photo- 
metric period of this phenomenon. In this work we refer to 
the fractional period offset of the common superhumps as 
the period excess e+, and the fractional period offset of the 
nodal/negative superhumps as the period deficit e_. 

In Paper II we presented visualizations and ray-traced 
simulation light curves which confirmed that the dominant 
physical source of the negative superhump signal is the 
changing depth of the accretion stream bright spot in the 
potential well of the white dwarf primary. The bright spot 
of an untilted disc will always be located on the outer edge 
of the disc, but if the disc is tilted the bright spot can im- 
pact the face of the disc. Consider that for a disc which is 
tilted and fixed in the inertial frame, the bright spot will 
impact the edge of the disc exactly twice per orbit, and be- 
tween these times will sweep across first one face of the disc 
and then the other. Assuming there are no dynamically- 
significant magnetic fields, the accretion stream will flow in 
the orbital plane of the binary, and impact the disc face or 
disc rim roughly along the line of nodes. Negative super- 
hump maximum occurs when the bright spot is at minimum 
radius - having fallen further in the primary star's poten- 
tial well, there is more specific kinetic energy to dissipate, 
and the bright spot is more luminous. The ray-traced light 



curves in Paper II show that negative superhumps observed 
from opposite sides of the disc will differ in phase by 180°. 
For a tilted disc that precesses slowly in the retrograde di- 
rection with a period PZ xec , the negative superhump signal 
has a period slightly less than the orbital period, 



1 

~pl 



Porl 



+ 



(2) 



2 NUMERICS AND PARAMETERS 

We simulate accretion disc dynamics us ing our 
smoo th ed particle hydrod y namic s (SPH) codj'] | Simpsonl 



19951; ISimpson fc Wood 



2007 



.19981 ; IWood fc BurkL _ 
Dolence. Wood, fc SilverF 20081 ). The work presented here 
was calculated in parallel with a large parametric study of 
common superhumps to be published separately. 

SPH is a Lagr angia n fluid dyn amics technique origi- 
nally developed by iLucvl (|l977l ) and iGingold fc Monaghanl 
(1 19771 ) that has s een wi de use in astro ph ysical 



applications (iBenzl Il990l; iMonaghanl |l992l. |2005|; 



2005; Wctzs tein et all 



Whitehouse. Bate, fc Monaghanl 
20081 ). The fluid properties are calculated by using a kernel 
interpolation in the scattered grid defined by the particles 
themsel ves. As is typical, we use a normalized spline 
kernel l|Monaghan fc Lattanziol Il985l ) that approximates 
a Gaussian but is identically zero beyond an interparticle 
distance of 2h, where h is the smoothing length. 

In their general forms, the momentum and en- 
ergy equations including gravitational bod y forces are 
ISimpson fc Woodll 19981 ; IWood fc Burke1l2007l ) : 

<fr 

dt 2 



VP 



+ fvis 



GMi 



-Tl 



GM 2 
r 2 3 



-r 2 , 



(3) 



da 
~dt 



= V ■ V + Cvisc, 

P 



(4) 



where f v i sc is the viscous force, e v i sc is the energy generation 
from viscous dissipation, n,2 = r — rin,M2 are the displace- 
ments from the stellar masses Mi and M2 , respectively, and 
the remaining symbols are defined as usual. 

The form of the SPH momentum equation for a particle 
i is given by 



d 2 r, 

dt 2 
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Ar + ^r I il + 1I.,)V,H'„ 



-Til 



GM 2 
ri2 3 



-r i2 , 



(5) 



where mj is the mass of particle j and W is the SPH kernel 
function. We u s e the artificial viscosity approximation of 
lLattanzio et al.l (|l986l ) evaluated between particles i and j, 



Ha = 



where 



-a/Mj + Vij ■ Tij < 0; 



otherwise; 



(6) 



1 A user-friendly version of this incorporating a graph- 
ical user interface and double-buffered graphics is avail- 
able for free download (www . astr o . f it . edu/wood/f itdisk.html; 
IWood. Dolence. fc Simpsonll2006l) . 
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where V; 



Vj 



(7) 

Ti — Yj. The the average of the 
sound speeds for particles i and j is c Sl y. For the simulations 
we present here, we adopt the parameters a = 1.5, /3 = 0.0 
or 0.5, rj = O.lh. As is typical in artificial viscosity prescrip- 
tions, only approaching particles feel a viscous force. We use 
a gamma-law ideal gas equation of state P — (7 — l)pu and 
use 7 = 1.01 since we do not include r adiative los ses. Within 
this model the sound speed is c s = */ 7(7 — l)u. 

In our simulations, all particles have the same smooth- 
ing length within a given simulation, varying smoothly from 
h = 0.0040a for the q = 0.75 run to h = 0.0067a for the 
q — 0.01 run, where a is the semimajor axis of the binary 
and a = 1 in system units. The SPH particles can have 
their own individual time steps as short as needed to re- 
solve the local dynamics. The largest timestep a particle can 
have is 5t° = P O rb/200, and other available timesteps are 
8t k = St°/2 k , where 1 < k < fc max and t ypica lly fc max = 9. 

As discussed in ISimpson fc Woodl (l998h . we integrate 
the internal energies using an action-reaction principle which 
is formally equivalent to the standard SPH energy equation. 
We assume that the sum of the changes in the internal en- 
ergies of all the particles over the previous time step (Sto) 
is directly proportional to the change in the bolometric lu- 
minosity of the disc over that same time interval, but do 
not include radiative transfer explicitly. With this assump- 
tion we can estimate a "simulation light curve" as the time 
series over timesteps n, 



L n = 



3 



fduY 
\dt). 



St , 



(8) 



and it is the frequencies present in these light curves that 
we discuss below. The dominant source of increase for the 
particles' internal energies is viscous dissipation, but we also 
include positive and negative P dV work in the summation. 

We begin our discs ab initio, injecting at a rate of typ- 
ically Pi n j = 2000 or 3000 particles per orbit until we reach 
our desired total of 100,000 particles in the disc. Our code 
has two options for handling particle accretion and injec- 
tion. The first option, and the one we use most often, is to 
promptly reinject at LI any particles that are ejected from 
the system or accreted onto Mi or M 2 . The second option 
decouples the injection rate from the accretion rate; injec- 
tion through LI is at a constant rate, and accreted or lost 
particles are not promptly replaced at the LI region. In the 
simulations we present in this work, we promptly re-inject 
during the initial evolution, so the total particle number 
simply increases linearly in time until the maximum num- 
ber of particles have been injected, and from that point on 
the particle number is is held constant. 

( Rt t< iVmax/P 

N=l (9) 
I A r max otherwise. 

Because the discs are initially forced to have a constant total 
number particles as a result of the prompt particle reinjec- 
tion after accretion or loss, the discs settle into a dynamical 
equilibrium. Discs that are not unstable to positive super- 
hump oscillations reach stationary states in the co-rotating 
frame. Those that are driven to superhump oscillations settle 
into a state of dynamical equilibrium where the dissipation 



light curve pulse shape is stationary from one superhump 
cycle to the next. 

To scale our models to physical units, we use t he empir- 
ical mass-radius relation of lPatterson et all (|2005h . The au- 
thors note that the data are "too sparse to yield trustworthy 
results" for secondaries with masses below O.O6M0. How- 
ever, needing an analytical relation in this regime whether 
trustworthy or not, we use a by-eye linear fit to the few data 
points shown in their Fig. 12 to obtain the modified mass- 
radius relation that we use to obtain the periods given in 
Table 1 below: 



O.23(M 2 /M ) 



M 2 < 0.06 M 



© 



R2 



O.62(Af 2 /M ) - 61 0.06 < M 2 /M < 0.20 (10) 



O.92(M 2 /M ) 



M 2 > 0.20M, 



©• 



Once the disc in a given simulation reaches a state of 
(dynamical) equilibrium, we impose a disc tilt of 5° on the 
particle positions and velocities and restart the simulation, 
evolving parallel sequences both with prompt replacement 
and with the mass injection rate set to zero. 



3 RESULTS 

3.1 Pre-tilt Evolution 

We simulated systems with mass ratios from q = 0.01 to 
q = 0.75. Simulations throughout this range yield negative 
superhumps in the light curves if the disc is tilted out of the 
plane and the simulation continues using prompt replace- 
ment of accreted/lost particles (see, e.g., Papers I & II). In 
Fig. \T\ we show "snapshots" of systems of 4 representative 
mass ratios at selected orbits and scaled to physical dimen- 
sions using the modifi ed secondary mass-radius relation of 
iPatterson et all (|2005l ) as given in Eq. (10) above. Our range 
of mass ratios includes systems with 0.03 < q < 0.35 that are 
unstable to positive superhump oscillations. In the Figure 
we show a disc of q = 0.05 that is near superhump maxi- 
mum, a disc of q — 0.20 that is near superhump minimum, 
as well as discs of q = 0.40 and 0.75 that are stable. The 
dimensionless disk thickness for these simulations is roughly 
H/r ss 0.04 as determined by-ey e using a plot of z versus 
r. For comparison, see Fig. 11 of I Wood et ail (|2005l ) which 
shows z vs. r for a simulation of q — 0.66 appropriate for 
that study of DQ Herculis. 

As an example of a simulation unstable to common su- 
perhump oscillations, we show in Fig.[5]the light curve of 200 
orbits of a q = 0.05 simulation for which the initial injection 
rate was 2000 particles per orbit. Common superhump oscil- 
lations begin at approximately orbit 70, and the disc settles 
into a state of dynamical equilibrium by roughly orbit 100. 
The average pulse shape and Fourier transform of orbits 100 
to 200 are shown in Figs. [3] and 2] respectively. The third 
harmonic dominates the average pulse shape, and this is re- 
flected in the Fourier transform. We will present the details 
of our apsidal superhump study in a separate publication. 

As an example of a system that is not unstable to ap- 
sidal superhump oscillations, we show in Fig. [5]the first 200 
orbits of a simulation with mass ratio q — 0.40. The initial 
particle injection rate for this simulation was 3000 per orbit, 
so 100,000 particles have been injected by orbit 34. By orbit 
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q = 0.05 (SH + max) 



0.20 (SH + min) 



0.40 (stable) 



0.75 (stable) 



Figure 1. Snapshots of mass ratios q = 0.05 (shown near su- 
perhump maximum) The centre of mass is indicated by a plus 
(+) sign. The primary white dwarf is shown to scale (J?wd = 
0.011 Rq for 0.8 Mq), and the Roche lobes are also shown. The 
snapshots are at orbital phase zero. The x axis points towards the 
bottom of the page, the y axis points to the right, and the z axis 
points out of the page. 
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Figure 2. The q = 0.05 simulation results from start to orbit 
200, by which time the disc is in a state of dynamical equilibrium 
defined as having common superhump oscillations that produce a 
light curve that repeats from one cycle to the next. In the figure 
above we sum 5 points in the light curve to yield 40 points per 
orbit, but no smoothing has been applied. Note that the scale of 
the lower panel differs from that of the upper panel by a factor 
of 2. 
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200, the disc has reached a stationary state in the co-rotating 
frame. Because we note some small trend in orbits 100-200, 
we ran the sequence out to orbit 400 and repeated some of 
the experiments. We found no difference in the resulting e_ 
or simulation light curves. 

The cause of the disc tilts in CVs is not com- 
pletely understood at the time of this writing. While 
in the low-mass X-ray binaries (LMXBs) it is believed 
that radiation pressu re can p r ovide the force neces- 
sary to tilt the disc (IPettersonl 19771: Iping fc Pettersonl 
19901: IPringldl 199rj. Il997l-| Foulkcs, Hasw ell. fc Murravll200ri 



Ivanov fc Papaloizou bOOa), this mechanism cannot explain 
the d isc tilt in CV discs. iBarrett. O'Donoghue. fc Warnerl 
(1988) suggested in their early study of TV Col that mag- 
netic fields on the secondary might be able to deflect the 
accretion flow through the LI region out of the orbital 
plane, thus breaking that symmetry, but as discussed in Pa- 
per II, as long as the deflection angle is stationary in the 
co-rotating frame and the mass transfer rate constant, the 
orbit-averaged angular momentum vector of the accretion 
stream is still parallel with the orbital angular momentum 
vector. Thus no tilt could be generated. iMurrav et al.l (|2002| ) 
showed that a disc tilt could be generated by instantaneously 
turning on a magnetic field on the secondary star of the sys- 
tem. This tilt decayed with time, but suggests that changing 
magnetic field geometries on the secondary could lead to disc 
tilts. The model which is most generally accepted is that the 



Figure 3. The average apsidal superhump pulse shape for orbits 
100 to 200 of the q = 0.05 simulation. We fold the simulation 
light curve on the fundamental period obtained using the Fourier 
transform. We bin into 50 samples per cycle, and show 2 cycles 
of the oscillation. Note that the third harmonic is a dominant 
component to the pulse shape. 
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Figure 4. The Fourier amplitude spectrum of orbits 100 to 200 of 
the q = 0.05 simulation. The peaks in the transform are labeled 
with their frequencies in units of cycles per orbit. As found in 
some observed systems, the fundamental (first harmonic) peak is 
not the highest-amplitude peak in the amplitude spectrum. 
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Figure 5. The q = 0.40 simulation results from start to orbit 
200, by which time the disc is in a stationary state. We sum 5 
points in the light curve, for a sampling of 40 per orbit. The scale 
of the lower panel differs from that of the upper panel by a factor 
of 2. 



primary white dwarf has a magnetic field, which can tilt the 
inner disc, and this information propagates outward in ra- 
dius, leading to a general disc tilt. We have completed some 
preliminary simulations which demonstrate that the non- 
parallel component (i.e., not parallel to the binary angular 
momentum vector) of the angular momentum of a tilted 
inner disc can propagate outward in radius and result in 
negative superhumps. 



3.2 Post-tilt Evolution 

We are working on simulations that explore the origin of 
the disc tilt in negatively superhumping CVs, but for this 
work, we simply impose a 5° tilt on a disc that has reached 
a state of dynamical equilibrium, and follow the evolution 
beyond that point. In Fig. [6] we show the side-view evolu- 
tion sampled every 5 orbits for sequences with mass ratios 
of q = 0.05, 0.40, and 0.75. The initial snapshot in each 
sequence shows the disc immediately after the disc tilt has 
been applied. Note that while the q = 0.40 and 0.75 models 
are outside the mass ratio range that yields positive super- 
humps, the q — 0.05 model shown is oscillating with positive 
superhumps as well as precessing in the retrograde direction. 

For use in discussions below, we define the direction 
"up" to be in the direction of the orbital angular momentum 
vector (it is the +z direction in our simulations). This also 
defines the terms "above the plane" and "below the plane." 
The geometry shown for the starting conditions at orbit 200 
in Fig. [5] has the line of centres of the stars perpendicular to 
the line of nodes, and the disc particles nearest the secondary 
are above the orbital plane. We define this to be negative 
superhump phase zero. At negative superhump phase 0.5, 
the line of centres and line of nodes are again perpendicular, 
but the particles nearest the secondary are below the plane. 
At negative superhump phases 0.25 (0.75) the line of centres 
and line of nodes are coincident, and the particles nearest the 
secondary have negative (positive) z velocities on average. 

For each simulation run, we typically continue the evo- 
lution of the titled disc for an additional 50 orbits for the 
two cases of (i) prompt particle replacement and constant 
total number of particles in the simulation, or (ii) no parti- 
cle replacement and a decreasing total number of particles in 
the simulation. We show sample simulation ( "bolometric" ) 
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Figure 6. The side-view evolution sampled every 3 orbits for 
simulations with mass ratios of q = 0.05, 0.40, and 0.75, from 
left to right. The Figure shows both that the nodal precession 
frequency is an increasing function of mass ratio and that the 
tilt is persistent. Note that the q = 0.05 simulation is oscillating 
with common superhumps as well as precessing in the retrograde 
direction. The plot limits in x of each subpanel are —0.8 to 0.8 
in system units where a = 1. The aspect ratio is 1:1 in the Fig- 
ure and only 2000 randomly-selected particles are shown for each 
snapshot. 



light curves for the post-tilt evolution of these two cases for 
the q — 0.40 run in Fig. [7] The upper curve of the Figure 
is the prompt-replacement case. The increase in mean lumi- 
nosity is the result of increased accretion onto the primary 
as a result of the accretion stream depositing new mass di- 
rectly to annuli of small radius. Because accreted particles 
are promptly replaced, this also enhances the particle injec- 
tion rate at LI over the pre-tilt average. The system relaxes 
to a new quasi-equilibrium after about 10 orbits. The lower 
curve of the Figure shows the case where we have shut off 
the injection of particles at LI completely, which we call 
the "no-stream" case. Here the system luminosity decreases 
both because there is no accretion stream bright spot, and 
because the total number of particles in the disc is decreas- 
ing monotonically. At orbit 450, only 73,162 of the original 
100,000 particles remain in the disc. The overall disc ge- 
ometry does not change significantly over t he 50 o rbits - 
e.g., the root-mean-squared radius r rms = \J — r~ of the 
particle distribution is r rms = 0.35a for both - but because 
each particle has fewer nearest neighbors with which it in- 
teracts, the total viscous dissipation as a function of time is 
decreased. 

In the inset of Fig.[7]we show the fractional amplitudes 
of the first 10 orbits of the post-tilt evolution for the two 
cases. To obtain the fractional amplitudes, we divide the 
original (200 samples per orbit) light curve by a copy of the 
light curve that has been boxcar-smoothed with a window 
width of 2 orbits and subtract 1.0. The smoothed curve will 
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Figure 7. The post-tilt light curves for the q = 0.40 simulation. 
We show both the light curve for the case of prompt replacement 
of accreted/lost particles (solid circles/upper curve) and for no 
replacement (open circles/lower curve). Here we have summed 10 
timesteps for 20 samples per orbit. To bring out the oscillations 
in the light curve, the inset figure shows the first 10 orbits of the 
same curves as fractional amplitudes as discussed in the text. The 
negative superhump signal is clearly present in the upper curve 
of the inset, but difficult to discern in the lower curve. 



Figure 8. The post-tilt light curves for the q = 0.05 simulation. 
We show both the light curve for the case of prompt replacement 
of accreted/lost particles (solid circles/upper curve) and for no 
replacement (open circles/lower curve). Here we have summed 10 
timesteps for 20 samples per orbit. The inset figure shows the first 
10 orbits of the same curves as fractional amplitudes as discussed 
in the text. The negative superhump signal beats with the positive 
superhump signal in the upper curve of the inset, but only the 
positive superhump signal is present in the lower curve. Positive 
superhumps continue long after the accretion stream is shut off, 
but comparison with Fig. [3] shows the pulse shape has changed. 



contain no trace of any strictly periodic signals that have an 
integer number of cycles per orbit, and will essentially free of 
the negative and positive superhumps signals. The residual 
fractional amplitude light curve will retain only frequencies 
higher than ~^ or b- Finally we then average groups of 10 
points to yield the plotted light curve which effectively has 
20 "exposures" per orbit. 

In the upper light curve of the inset of Fig. [7J the nega- 
tive superhump signal can be clearly seen, while in the lower 
curve it is not obvious above the noise. Note that the signal 
in the upper curve has a frequency of approximately 2 cy- 
cles per orbit. As discussed in Paper II, the reason for this is 
that the dominant source of the negative superhump signal 
is the transit of the accretion stream bright spot across the 
face of the tilted disc. Because our bolometric light curves 
reflect the viscous dissipation of all the particles in the sim- 
ulation, there are negative superhump maxima correspond- 
ing to the bright spot sweeping across both faces of the disc 
in the prompt-replacement light curve. The minima in this 
light curve correspond to the times at which the accretion 
stream impacts the disc rim along the line of nodes. In an 
actual system, an observer only sees the effects of the bright 
spot sweeping across a single face of an optically-thick disc, 
and so sees a single maximum and single minimum per orbit. 

For comparison, we show in Fig. |H]the post-tilt evolu- 
tion of the q — 0.05 simulation, which has developed positive 
superhumps before the tilt is applied at orbit 200. In this 
figure, we offset the no-stream light curve downward for clar- 
ity. Note that the upper curve (with-stream) shows evidence 
for beating between the negative and positive superhump 
signals, whereas the lower curve (no-stream) shows only the 
positive superhump signal with no evidence for beating with 
a negative superhump signal. Again, the overall luminosity 
declines because the total particle number and hence mean 
number of nearest neighbors per particle in the simulation 
is decreasing, which decreases the total viscous dissipation. 
We note, however, that the positive superhump amplitude 
is nearly constant over the 50 orbits even though there is no 
accretion stream and hence no accretion stream bright spot. 



Thus, the primary contribution to the positive superhump 
signal is viscous dissipation within the disc itself, and not 
the varying depth around the orbit of the accretion stream 
bright spot in the potential of the primary white dwarf as is 
commonly (mis)stated in the literature. 

In the inset of Fig. [8] we show the first 10 orbits of 
the two cases as fraction amplitudes, calculated as detailed 
above. In the upper light curve, there are clearly two signals 
present as the light curve pulse shape evolves during the 
10 orbits. In the lower curve, only the positive superhump 
signal is present, and the fractional amplitude appears to 
increase during these 10 orbits. A comparison with the q — 
0.05 pre-tilt average pulse shape (Fig. [3} does show that the 
harmonic complexity of the light curve has decreased as the 
power in the third harmonic is noticeably reduced. 

3.3 Amplitude Spectra and Pulse Shapes 

A primary goal of this work is to give analytical expressions 
for the negative superhump period deficit e_ as a function 
of mass ratio q and orbital period P or b- To achieve this, 
we simulated systems with mass ratios spanning the range 
0.01 < q < 0.75. We evolved the systems to a state of equilib- 
rium, tilted the discs by 5° with respect to the orbital plane, 
and continued the simulations either by promptly replacing 
the accreted/lost particles, which we refer to as the "with 
stream" (ws) case, or by shutting off particle injection at LI 
completely, which we refer to as the "no stream" (ns) case. 
Simulations were continued for typically 50 orbits. For each 
simulation, we take a Fourier transform of the simulation 
bolometric light curve converted to fractional amplitudes as 
described above, and after removing spurious points that 
are more than 6cr away from the mean. After calculating 
the Fourier transform, we interpolate to find the frequen- 
cies and amplitudes of the 50 highest peaks. The peaks in 
the transform have a full width at half-maximum of roughly 
8v ~ 1/50 cyc/orb. 
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In Fig. [9] we show for representative mass ratios the 
Fourier transforms of the post-tilt bolometric light curves 
in both the with-stream and no-stream cases. The negative 
superhump signal in the light curves is reflected in the ampli- 
tude spectrum as a peak at a frequency slightly greater than 
2 eye orb -1 . This peak can be found in all spectra where the 
accretion stream is present, and in the spectra with mass 
ratios of q — 0.30 and above in the no-stream case. The 
q — 0.05 post-tilt spectra show differences from the pre-tilt 
spectra shown in Fig. [4] above, reflecting the differing pulse 
shapes. Note that the post-tilt positive superhump first- 
harmonic (fundamental) amplitude is significantly higher in 
the no-stream case as compared with the with-stream case. 
This result is even more dramatically shown in the q = 0.20 
amplitude spectra. The positive superhump amplitudes can 
increase in the no-stream case because the disc is no longer 
accreting mass of low specific angular momentum at the 
disc edge, and so can expand slightly and be driven more 
strongly. Careful comparison with Fig. [4] shows that ampli- 
tude of the the first harmonic is also stronger in the with- 
stream case of the post-tilt evolution. 

In the q — 0.20 with-stream amplitude spectrum, there 
is power at a frequency just above 1.0 eye orb -1 , but this 
power does not indicate a physical oscillation of the disc. 
Rather, this is a combination frequency between the first 
harmonic of the positive superhump frequency and the neg- 
ative superhump signal, resulting from non-sinusoidal pulse 
shapes for the two signals. Inspection of the peak frequencies 
in the figure shows that the frequency difference A between 
the simulation negative superhump signal at 2za_ and the 
second harmonic of the positive superhump signal at 2v+ 
yields combination frequencies of nv+ + A for several har- 
monics of the positive superhump fundamental. 

Note that the amplitude spectra of simulations lack- 
ing positive superhumps have no significant power near 
v = 1 eye orb -1 because of the inherent symmetry discussed 
above. Specifically, the transiting of the bright spot releases 
equal amounts of power on both faces of the disc, and the 
same must also be true of the tidal effects that give rise 
to the component of the negative superhump signal that 
arises even in the absence of an accretion stream. This can 
also be demonstrated by plotting the average pulse shapes 
for negative superhumps, which we show in Fig. QIJ] for a 
more limited but still representative sample. To compute 
the average pulse shapes, we compute the fractional ampli- 
tude residuals as described above, fold the light curves on 
P- — Iff- found from the Fourier transform, and bin into 
40 bins per cycle. Because we subtract the time of the first 
data point from the time series before folding, zero phase 
in these curves is as defined geometrically above (with the 
line of centres perpendicular to line of nodes, and particles 
nearest secondary above the orbital plane). For the mass 
ratios shown, it is clear that to within the noise, the sim- 
ulation bolometric pulse shape rep eats twice per nega tive 
superhump cycle. Note finally that I Montgomery! (|2009| ) in 
her study of negative superhump simulations claims that the 
half-cycle average pulse shapes differ (her Figs. 8 and 13). 
However, her pulse shapes are undersampled and appear to 
contain an odd number of bins per cycle. Her Fourier trans- 
form shows no significant power near v = 1, which would 
have to be present for her claim to be valid. 

The character of the bolometric average pulse shape 
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Figure 9. Fourier transforms of selected post-tilt simulations. In 
the figure we show Fourier transforms of 50 orbits of simulations 
with q = 0.05, 0.20, 0.40, 0.55, and 0.75. After imposing the disc 
tilt, we continued the simulations either with prompt replacement 
of accreted/lost particles (i.e., "with stream," abbreviated "ws"), 
or with mass injection shut off completely (i.e., "no stream", ab- 
breviated "ns"). The Fourier transforms shown are all on the 
same scale, but the curves have been offset for clarity. See text 
for discussion. 



changes as a function of mass ratio. In the no-stream case, 
the amplitude decreases significantly from the q = 0.40 sim- 
ulation to the q = 0.55 simulation, and is consistent with 
zero to the eye for the q = 0.75 simulation (not shown). In 
the with-stream case, the q = 0.40 average pulse shape has 
two rounded maxima per cycle, whereas the 0.75 average 
pulse shape is saw-toothed. The difference may lie in the 
fact that the accretion stream has less specific angular mo- 
mentum for increasing q, and more specific kinetic energy 
at the bright spot (see again Fig. [IJ. 

We present the results of the Fourier analysis in Table 1, 
where we give for each run (defined by q and viscosity coeffi- 
cients a and 0) the pre-tilt positive superhump frequency h>+ 
(in units of cycles per orbit) and the orbits for which that 
was calculated, and the post-tilt orbit range and negative 
and positive superhump frequencies and fractional ampli- 
tudes for both the with-stream and no-stream cases. Entries 
marked with a dash (-) indicate no significant periodic sig- 
nal was found. We note that the tabulated values for the 
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Figure 10. Average pulse shapes for mass ratios q = 0.40, 0.55, 
and 0.75, computed from the bolometric simulation light curves 
and for both the with-stream and no-stream cases (the q = 0.75 
no-stream average light curve has negligible amplitude and is not 
shown) . The curves show that the half-cycle pulse shapes are iden- 
tical within the noise. The character of the pulse shapes changes 
as a function of mass ratio, becoming more triangular for higher 
mass ratios for the with-stream case. 



negative superhump frequencies are half that found in the 
Fourier transforms of the bolometric light curves. The frac- 
tional amplitudes are those of the peaks 2v- peaks present 
in the transforms (see, e.g., Fig. [9}. 

There are a few items of note from this table. First, ev- 
ery model with steady accretion yields a negative superhump 
signal if the disc is tilted, and the higher the mass ratio, the 
higher the negative superhump frequency in cycles per orbit. 
The fractional amplitudes of the negative superhumps in the 
with-stream case varies within a range of roughly ~0.6 to 
1.4%, and those in the no-stream case vary within the range 
~0.17 to 1.1%. The fractional amplitudes in the with-stream 
case do not appear to vary in a well-defined way in the range 
of q where positive superhumps are present simultaneously, 
but they grow monatonically with q for simulations without 
simultaneous superhumps, as already discussed in reference 
to Fig. [10] 

In the absence of an accretion stream, only systems with 
mass ratios q > 0.30 display a negative superhump signal. 
The fractional amplitude of the no-stream negative super- 
hump signal rises rapidly from zero for q<0.25 to ~1% for 
q = 0.35, and then falls more gradually to q = 0.75, at which 
point the peak in the amplitude spectrum is barely above 
the noise. We propose that the physical reason for this is 
a forced driving of the tilt instability mode initially recog- 
nized bv lLubowl (|1992| ). Near q = 0.30, most of the mass in 
the simulated discs is near the 3:1 resonance orbit, and even 
at q = 0.40, the outer disc particles complete ~3 orbits per 
negative superhump cycle, thus periodic disipation can take 
place. This driving becomes increasingly weaker as the outer 
disc particles move further from the resonance orbit. Note 
that although for q ~ 0.35 the fractional amplitude of the 
no-stream case can be larger than that of the with-stream 
case (which must of course also include the no-stream con- 
tribution), this is simply the result of division by the lower 
mean light levels for the smoothed no-stream light curves 
(see, e.g., Fig.0. 

Systems in the mass-ratio range 0.03 < q < 0.35 are un- 
stable to positive superhump oscillations. We have run a 
q = 0.025 simulation out to orbit 1000 and no superhumps 



developed, confirming t he lower bound we determined in 
ISimpson fc Woodl (|1998l ) . Interestingly, but not surprisingly, 
the upper bound appears to be a function of mass ratio. 
The tabulated results show that for q = 0.35, the a — 1.5, 
(5 = 0.0 simulation, positive superhumps develop, whereas 
they do not develop for the higher-viscosity a — 1.5, ft = 0.5 
simulation which was run out to orbit 500. However, once 
the equilibrium disc at orbit 200 is tilted, we find positive su- 
perhumps clearly developing in the no-stream post-tilt evo- 
lution at roughly orbit 240, resulting in a clear peak in the 
Fourier transform at v+ ~ 0.82 cycles per orbit. The with- 
stream result is less convincing. There is a peak at the same 
location in the amplitude spectrum, but it is only about 
twice the height of the surrounding noise. Finally, we note 
that in most cases the results for the frequencies of the nega- 
tive or positive superhumps are not strongly affected by the 
choice of viscosity over the range we have used. 

3.4 Ray- Traced Visualizations 

In order to understand better the source of the nega- 
tive superhump modulation, we visualized and animated 
our results using IDL object graphics. We show still 
frames from these animations in Fig. 1111 These and 
other accretion disc animations can be viewed online 
at www.astro.fit.edu/wood/visualizations.html. In the 
figure, we include representations of the Roche-lobe profiles 
in the orbital plane, the primary star approximately to scale, 
the line of nodes, and color-code the particles by the lumi- 
nosity over the previous timestep. The light curves shown 
are ray-traced as described in Paper II - we sum the energy 
generation (Eq. [5} over only the "surface" particles along a 
given line of sight. The top light curve in the figure corre- 
sponds to the i = 0° light curve (i.e., as viewed from the 
positive z axis), and the bottom curve corresponds to the 
i — 180° light curve. Note that these are anti-phased with re- 
spect to each other for both the with-stream and no-stream 
cases, and that the ray-traced light curves only show 1 max- 
imum per superhump cycle. The change in color of the light 
curves indicates the current time for that frame. 

The images in the left column of the figure result from 
the simulation that includes the stream, and shows 1 com- 
plete orbit. The right column shows the same orbit for the 
simulation without the stream. Note that the time value 
shown in the upper right hand corner of each frame is the 
time in orbits, not in superhump phase. The impact region 
of the accretion stream is visible in the with-stream images 
at times 402.45 and 403.45 (top and bottom panels on the 
left hand side). 

3.5 Fits to the Tabulated Results 

In Fig. Q2] we show the negative superhump period deficit 
and positive superhump period excess as functions of the 
mass ratio q. The Figure also includes fits to the numerical 
results. In fitting the negative superhump results, we find 
that the following four-parameter fit provides an excellent 
approximation to the numerical results over the range 0.01 < 
q < 0.75. The positive superhumps are not the focus of this 
publication, but we include a preliminary fit to these results 
for completeness: 
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a Calculated using our adapted version of the empirical mass-radius relation of lPatterson et all ll2005f ) given by Eq. dlOI I 
above. 

" Frequencies are given in units of cycles per orbit (cyc/orb). 
c Blended peak. 

^ Signal grows to observable amplitude at roughly orbit 240. 



e_ = -0.02263<j 1/2 - 0.277q + 0.471<? 3/2 - 0.249? 2 (11) 

e+ = 0.238g + 0.357? 2 (12) 

Because it is more common that q is the desired quan- 
tity for a system with observed superhump and orbital peri- 
ods, we also give fits to q versus |e_| and e+. We note that 
the apsidal precession rate depends somewhat on the disk 
sound speed, and that this effect becomes more important at 
smaller q. We will explore this more fully in our forthcoming 
parametric study of apsidal superhumps. 

q = -0.192|£_| 1/2 + 10.37|e_|-99.83|e_| 3/2 +451.1|e_| 2 (13) 

q = 3.733e+ - 7.898e+ (14) 

The superhump period excess for the q = 0.10 simu- 
lation appears to deviate substantially from the trend sug- 
gested by the remaining points, and we originally suspected 
that there was a problem. However, we have checked inde- 
pendent simulation runs with differing values of the viscosity 
parameters (e.g., a = /3 = 0.5), smoothing length, and total 
particle number, and they consistently group near the plot- 
ted point. We will discuss possible reasons for this in our 
forthcoming publication presenting the results of our pos- 
itive superhump parametric study. This is one reason why 
we emphasize that Eq. Q2] is preliminary. 

In Fig. [T3] we plot the magnitude of the ratio (j> = 
|e_/e+| of the negative superhump period deficit to the 
positive superhump period excess over the mass ratio range 
0.03 < q < 0.35. As expected from Fig. 1121 the ratio is near 
unity for the smallest mass ratios, and declines as mass ratio 
increases, a result that was hinted at in Fig. 2 of lRetter et al.l 



(|2002i b Again in Fig.[T3]we find that the q = 0.10 point de- 
viates from the trend indicated by the remaining points. 
Because of this we, offer 3 different linear fits to the results. 
Fitting all the data points we obtain 

<t> = 1.033 - 2.16g. (15) 

If we chose not to include the q = 0.10 point, then the linear 
fit is 

4> = 1.091 - 2.34<7. (16) 

Finally, in consideration of the fact that negative super- 
humps are primarily only observed in longer-period systems 
(i.e., above the period gap), we have a fit to the points in 
the range 0.10 < q < 0.35, 

<t> = 0.837 - 1.37<?. (17) 

Up to this point we have presented our results as a 
function of mass ratio, because that is the fundamental pa- 
rameter in the simulations. However, the mass ratio is much 
more difficult to measure observationally than orbital or su- 
perhump periods, and so in Fig. [14] we show the negative 
(positive) superhump period deficit (excess) as a function 
of -Porb , where we used the modified iPatterson et al.l (|2005T 1 
secondary mass-radius relation (Eq. 10) to scale our mass 
ratios to periods, assuming that the white dwarf primary 
mass is M\ = 0.8 A/0. This empirical mass-radius relation 
yields a minimum P OI b, and for this reason we confine our 
fits to mass ratios of q > 0.10. These relations should help 
guide observers to determine if a putative photometric pe- 
riod detection is plausibly the result the negative super- 
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Figure 12. The superhump period deficit (solid circles) and 
excess (open circles) as a function of mass ratio q. The four- 
parameter fit to the negative superhump results is shown as a 
solid line, and the two-parameter fit to the positive superhump 
results is shown as a dotted line. Note that the points at q = 0.03 
and 0.05 are nearly degenerate (e_ as e+) in the Table and over- 
lap substantially in the figure. 
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Figure 13. The ratio of the superhump period deficit to excess 
as a function of mass ratio q. Because e+ for the the q = 0.10 
simulation appears to be an outlier, three different linear fits are 
shown, as discussed in the text. 




Figure 11. Panel figure showing one complete orbit each of 
the q = 0.40 simulations with (left panel) and without (right 
sequence) an accretion stream. The system is shown in the re- 
duced mass frame (fixed primary location). In the Figure, the 
intersection of the Roche lobes with the orbital plane is shown 
as a yellow line, the particles are color-coded by luminosity over 
the previous timestep, with the brightest particles rendered as 
slightly-larger blue points. The ray-traced simulation light curve 
calculated for i = 0° is shown above the disc in each frame, and 
below for i = 180°. Note these are antiphased. The portion of 
the lightcurve preceding the current timestep is a lighter color 
than the remainder. The line of nodes is shown as a blue bar 
in the orbital plane, and it precesses in the retrograde direction. 
These still frames are from animations that may be viewed as at 
www. astro . f it . edu/wood/visualizations . html. 



humps. Fitting all the negative superhump results over the 
range 0.10 < q < 0.75, we obtain 



e_ = -0.0142 - 0.00852P orb . 



(18) 



Because the negative superhump results appear nearly lin- 
ear over the range 0.10 < q < 0.25, we fit these and the 
remaining points separately 

r 0.01225 - 0.0237P O rb 0.10 < q < 0.25, 
e- = <^ (19) 
t -0.00722 - 0.0101P orb 0.30 < q < 0.75. 

Finally, the best linear fit to the positive superhump period 
excess over the range 0.10 < q < 0.35 is 

e+ = -0.0226 + 0.0415P orb . (20) 



4 OBSERVATIONAL RESULTS 

4.1 Systems with Main Sequence Secondaries 

In Table 2 we list 21 systems that have published negative 
superhump period deficits. For each system, we list the or- 
bital period, the published period deficit, and the published 
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E_ (0.10 < q < 0.75) 

e_ (0.10 < q < 0.25) 

E_ (0.30 < q < 0.75) 

(0.10 < q < 0.35) 




Table 2. Observational Results. 



Figure 14. The superhump period deficit (solid circles) and 
excess (open circles) versus orbital period calculated assuming 
Mi = 0.8Mq and the secondary mass-radius relation given by 
Eq. mot . The linear fits are discussed in the text. 
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Figure 15. Observed superhump period deficit (solid circles) and 
excess (open circles) versus orbital period for systems listed in 
Table 2. The fits are as presented in Fig. 1141 The scatter of the 
observations is substantial. We do not include the ultra-compact 
systems AM CVn or V1405 Aql in the figure as the secondaries in 
these systems are not main sequence stars, and so these systems 
follow a different q vs. P rb relation. 



period excess if it exists. The data are presented graphically 
in Fig. 1151 In the figure, the solid circles represent the neg- 
ative superhump period deficits, and the open circles the 
period excesses. While most of the reported period deficits 
are found near the best-fit relationship we find from our 
numerical results, a handful of the observed period deficits 
are quite discrepant. In general the systems require very long 
photometric time series observations over a period of a week 
or more and at multiple longitudes for best results. It may 
be useful to revisit several of these systems to determine if 
the published period deficits can be confirmed. 

Negative superhumps have been tentatively identified 
in a handful of additional systems, but w e did not include 
these systems in Table 2 and Fig. 1151 iTramposch et all 
|2005l ) suggest that SDSS J210014.12+004446.0 shows neg- 
ative and positive superhumps with periods of 1.96 ± 0.02 
and 2.099 ± 0.002, respectively. The scatter in their phase- 
resolved spectroscopy suggests an orbita l period near 2 hr, 
but the value is not well determined. Retter. fc Lul 
l|2005al ) suggested that V1193 Ori displayed negative su- 
perhumps in their photometric time series observations, 
but to conclude that they had to assume that the ac- 
tual orbital period was aid -1 alias of the published 
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-Porb(h) 


£_ 


£ + 


Refs 


AM CVn 


0.286 


-1.7 


2.2 


1, 2, 3 


V1405 Aql 


0.834 


-0.71 


0.90 


4, 5, 6 


V1159 Ori 


1.492 


-7.6 


3.2 


7, 8 


ER UMa 


1.528 


-7.5 


2.1 


9, 10 


IR GEM 


1.642 


-3.1 


5.2 


11 


V503 Cyg 


1.865 


-2.6 


4.3 


12 


BF Ara 


2.020 


-2.4 


4.5 


13 


RX J1643+3402 


2.893 


-3.0 




14 


AH Men 


2.950 


-2.4 


3.1 


15, 16 


V442 Oph 


2.984 


-2.8 
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TT Ari 


3.301 


-3.4 


8.5 


17, 18, 19 


V603 Aql 


3.314 


-3.0 


5.6 


20 


RR Cha 


3.362 


-2.7 


3.1 


21 


V751 Cyg 


3.467 


-3.5 




22, 23 


PX And 


3.512 


-3.0 


9.0 


3, 24 


V2574 Oph 


3.546 


-4.1 




25 


HS 1813+6122 


3.55 


-4.5 




26 


BH Lyn 


3.741 


-7.0 




27 


KR Aur 


3.907 


-3.5 




28 


SDSS J0407-0644 


4.084 


-2.4 




29 


TV Col 


5.486 


-5.5 


15 


30, 31 



Reference s llSkillman et alj (Il999l'>; 2|Pattersonl lll998h ; 3 
Patterson! Jl999h: 4lChou. Grindlav. fc Bloserl feOOll l: 5 
Hu. Chou, fc Chung! l|2008l); 6lRetter et alj <20o3) ; 7 



Patterson et al.l Jl995T > ; 8 iThorstensen et alj dl997h ; 9 lGao et all 
jl99Sh; lOlZhao et alj ll2006h ; ll lFu et al.l J2004l ; 12 
Harvey et al.l lll995h 



13 



Olech. Rutkowski. fc Schwarzenbere-Czernvl 
Patterson et al.l j2002T); 15|Pattersonl dl995T ); 16 



14 



Rodriguez-Gil, Schmidtobrcick, Gansickc (2007aj); 17 
Skillman et all dl998T); 18lAndronov et alj Jl999h; 19lWu et alj 
j2002h; 2Q|Patterson et alj [1199711; 2llWoudt fc Warner! (I 2002T) ; 
22 |Patterson et alj d200lT ); 23 IPapadaki et alj <200S|) ; 24 



Stanishcv et al 



-] J2002fl: 25lKang et al.l ||2006|) ; 26 
Rodriguez-Gil et alj [|2007bh; 27 " 
Istanishev, Kraicheya, fc Genkovl [|2006h ; 2 8 iKozhevnikovl ll2007l) ; 
29 lAk et alj [|2005bh ; 30 lHellierl dl993h ; 31 iRetter et al.l d2003h ; 



value of -Porb = 0.165 d (iRingwald. Thorstensen. fc Hamwevl 
1 19941 ; iPapadaki et~ai1 120041 ). IPapadaki et alj (|2006l ) were 
able to confirm the published value of the orbital period 
and were not able t o confi rm any of the periods reported by 
lAk, Retter. fc Liul l|2005al ). The long-period system TX Col 
was claimed to potentially show bo t h pos itive and negative 
superhumps bv lRetter. Liu, fc Bos! l|2005h based on single- 
site data, however they warn that their observations need 
confirmation. 



4.2 V1159 Ori and ER UMa 



IPatterson et al.l l|l995h show evidence of a 82.7-min period- 
icity in V1159 Ori which occured on two separate nights, 
and they identify this as a negative superhump signal. Us- 
ing the spectroscopically -determined P or b = 89.5363 min 
l|Thorstensen et akl [1997) , this yields a period deficit of 
— 7.6%, where we'd expect a value closer to e _ ~ — 1.6% 
based on our r esults. Similarly, iGao et all (| 19991 ) and 
IZhao et alj l|2006j ) report finding negative superhumps in ER 
UMa also with a very large period deficit of e_ = —7.5%. 
Here the period is found to vary night-to-night, and the run 
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lengths are short. If these signals can be confirmed, then 
they are perhaps the first observed signals at u+ + A dis- 
cussed above in reference to the q = 0.20 with-stream am- 
plitude spectrum (Fig. [9} . However, the absence for either 
system of a detected signal near what should be the "true" 
2u- weakens this interpretation. Additional observations are 
clearly needed. 



4.3 AM CVn & V1405 Aql 

AM CVn is perhaps the best studied of the hel ium-mass- 
trans fer binaries, and the prototype for the class l|Nelemansl 
120051 ) . The most extensive photometric study to date for this 
object is that of Skillma n et al.l (| 19991 ). The positive super- 
hump excess is e+ — 2.2%, which using Eq. (11411 above yi elds 
a mass ratio estimate of q ~ 0.08. ISkillman et al.l (1 19991 ) re- 
port e- = -1.7% for AM CVn, which from Eq. above 
yields an estimate of q m 0.06. These estimates are consistent 
with each other and with previous e stimates based the ob - 
served q versus e+ relati onship (e.g. iPatterson et al.lll995l ). 
However, more recently. iRoelofs et al.l 1 20061 ) presented the 
results of an extensive time-series spectroscopic study of AM 
CVn in which they find a kinematic feature that appears 
analogous to the "central spike" feature observed in some 
longer period AM CVn stars. If, as in these other systems, 
the observed central spike of the He II 4471 line traces the 
projected velocity and phase of the accreting white dwarf, 
then q = 0.18±0.01 for AM CVn (jRoelofs et al.ll2006D . which 
is substantiall y higher t han t he q inferred from the super- 
hump signals. IPearso 

3 <|2007h explores the possibility that 
helium discs in AM CVn systems may have pressure effects 
operating in their discs that differ significantly from those 
characteristic of hydrogen-dominated discs, and this cou- 
pled with possible differing se condary star formation chan- 
nels for these systems (see, e.g.. lRoelofs. Nelemans. fc Grootl 
l2007l ; lYunge"isonll2008l ) may preclude the existance of a well- 
defined superhump period excess to mass ratio relationship. 
However, other AM CVn stars with measured superhump 
and orbital periods generally appear to be consistent with 
the hydrogen-disc results, so additional effort is required to 
resolve these issues. 

V1405 Aql is a LMXB first discovered as an X-ray 
source exhib iting type-I burst s, suggesting a neutron-star 
primary fsee lRetter et al.ll2002l . and references therein). Ob- 
servations indicate an optical period about 1% longer than 
the 0.834-h X-ra y orbital per i od, co nsistent with common 
superhumps, and lRetter e t al. ( 2002]) report the detection of 
a periodicity 0.828 h, which they claim to result from neg- 
ative superhumps. Thus, for V1405 Aql the negative super- 
hump period deficit is e_ = 0.71% which implies q — 0.021 
using Eq. (I13|l above, and the positive superhump period 
excess is e+ — 0.90%, for which Eq. (|14p implies q = 0.033. 
For a 1.4 Mq neutron star primary, these results suggest 
a secondary mass in the range 0.03-0.05 Mq, and because 
we find positive superhumps only for q > 0.03, we would 
suggest the secondary star has a mass near 0.05 Mq. Given 
the discrepency for AM CVn between the spectroscopically- 
determined mass ratio and that inferred from the superhump 
periods, the above determinations must of course be taken 
as preliminary. 



5 CONCLUSIONS 

We report on the results of an extensive parametric study of 
negative superhumps using the method of smoothed parti- 
cle hydrodynamics. Negative superhumps are visible when a 
cataclysmic variable system has a disc which is tilted out of 
the orbital plane. Torques operating on the tilted accretion 
disc cause it to precess in the retrograde direction, at a rate 
that increases with system mass ratio q. 

Our results reveal that the negative superhump signal 
is dominated by the transit of the accretion stream impact 
spot across the face of the tilted disc. The deeper in the 
primary's potential well the impact point on the face of the 
tilted disc, the brighter the instantaneous luminosity. Simu- 
lations with mass ratios spanning the range 0.01 < q < 0.75 
and non-zero accretion rates through LI produce a negative 
superhump signal in the light curve. The signals are an- 
tiphased for observers on opposite sides of the orbital plane. 

The most surprising finding of this project was the dis- 
covery of the existence of a negative superhump signal in 
simulations for which the accretion stream had been shut off 
completely. This previously unknown and unsuspected effect 
is only present in simulations with q > 0.30. The signal has 
a fractional amplitude that peaks at ~1% for q ~ 0.35, and 
falls gradually with increasing q. We propose that the source 
of this signal is dissipation associated w ith forced driv ing of 
the 3 : 1 vertical resonance identified bv lLubowl (|1992| ). The 
driving is strongest for lower mass ratios as these have most 
of the mass (particles) in the region of the 3:1 resonance. 
The driving gets weaker with increasing q as the driving 
frequency pulls away from the response frequency. 

We give empirical fits to our simulation results which 
should be useful to observers working to identify negative 
superhumps in nature. While most published negative su- 
perhump period deficits are in reasonable agreement with 
our numerical simulation results, a handful are in serious 
disagreement. Because the physics of retrograde disc preces- 
sion as it applies to these systems is likely to yield intrin- 
sically small scatter about a well-defined trend, we suggest 
that it may be useful to observe the outlier systems with 
global campaigns in an effort to reduce the observational 
scatter. If, however, new observations unambiguously con- 
firm the large negative superhump period deficits reported 
for these systems, then perhaps the physics of photometric 
signal generation from tilted accretion discs in cataclysmic 
variable systems is even more interesting and complex than 
we currently know. 
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